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Abstract Solar flares observed in the 200-400 GHz radio domain may exhibit a 
slowly varying and time-extended component which follows a short (few minutes) 
impulsive phase and which lasts for a few tens of minutes to more than one hour. 
The few examples discussed in the literature indicate that such long-lasting 
submillimeter emission is most likely thermal bremsstrahlung. We present a 
detailed analysis of the time-extended phase of the 2003 October 27 (M6.7) 
flare, combining 1-345 GHz total-flux radio measurements with X-ray, EUV, 
and Ha observations. We find that the time-extended radio emission is, as 
expected, radiated by thermal bremsstrahlung. Up to 230 GHz, it is entirely 
produced in the corona by hot and cool materials at 7-16 MK and 1-3 MK, 
respectively. At 345 GHz, there is an additional contribution from chromo- 
spheric material at a few 10** K. These results, which may also apply to other 
niillimeter-submillimeter radio events, are not consistent with the expectations 
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from standard semi-empirical models of the chromosphere and transition region 
during flares, which predict observable radio emission from the chromosphere at 
all frequencies where the corona is transparent. 

Keywords: Radio Bursts, Association with Flares; Radio Bursts, Microwave; 
X-Ray Bursts, Association with Flares; Flares, Relation to Magnetic Field; 
Chromosphere, Active 



1. Introduction 



Since 2000, new instrumentation has allowed us to observe solar flares at sub- 
millimeter wavelengths. Such observations are routinely carried out by the Solar 



Submillimeter Telescope (SST; Kaufmann et al., 2008) at 212 and 405 GHz, and 
studies of ten major bursts detected by this instrument have been reported in 
the literature. Measurements at 230 and 345 GHz have been obtained for two 
more flares by the Koln Observatory for Submillimeter and Millimeter Astron- 



omy (KOSMA; Liithi, Magun, and Miller, |2004[ Liithi, Liidi, and Magun, |2004[ ) 
telescope, one of them being also observed at 210 GHz by the Bernese Multibeam 
Radiometer for KOSMA (BEMRAK). 

Above 200 GHz, these events exhibit an impulsive phase lasting for a few 
minutes, which is sometimes followed by slowly varying and time-extended emis- 
sion (the "gradual" phase). While, for some events, the >200 GHz emission 
appears as the extension toward high frequencies of the gyrosynchrotron emission 



seen in the microwave domain (Trottet et al., 2002 Liithi, Magun, and Miller, 



2004 Raulin et al., 2004 Gimenez de Castro et al, 2009|), other events exhibit an 

Such 



20041 



unexpected upturn towards the THz domain (e.g., Kaufmann et al.. 
spectra with positive slopes in the millimeter-submillimeter domain have been 
measured during both the impulsive phase (Kaufmann et al., 2004 Silva et al.. 



2007 Kaufmann et al, 2009) and the gradual phase (Trottet et al, 
Magun, and Miller, [2004jTuthi, Ludi, and Magun, [2004 ) 



2002 Liithi, 



Although various theoretical ideas have appeared in the literature (e.g., Kauf- 



mann and Raulin, 2006 Fleishman and Kontar, 2010), the emission mechanism 



responsible for a positive slope in the sub- THz domain remains uncertain, at 
least during the impulsive phase. During the gradual phase, it has been shown 
that thermal bremsstrahlung from the chromosphere and the corona may ac- 
count for the observed radio spectrum, particularly during the late decay of 
the flare. This is the case for the 2000 March 22 (XI. 1) flare (Trottet et al.. 



2002[ ) and for the 2003 October 28 (X17.2) flare (Trottet et a/., |2008| ), where the 
radio spectrum, late in the gradual phase, is roughly flat below 200 GHz and 
increases at higher frequencies. Such spectra were interpreted as optically-thin 
thermal bremsstrahlung from the corona below ~200 GHz, and it was suggested 
that the excess emission measured at higher frequencies arose from thermal 
bremsstrahlung from the lower atmosphere. 

The best example of a thermal phase at 230 and 345 GHz reported so far 
is the time-extended component of the 2001 April 12 (X2.1) flare, which lasts 
for more than two hours and during which the radio spectrum is flat first up to 
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230 GHz and later up to 89 GHz, with a positive slope above these frequencies 
(Liithi, Magun, and Miller, 2004 1. However, although the radio emission during 
the gradual phase has been conjectured to arise from separate contributions 
of coronal and chromospheric sources, the few available observations have not 
allowed a characterization of the plasmas involved. 

This paper presents an analysis of radio observations of the 2003 October 27 
(M6.7) flare at -12:30 UT obtained in the 1-405 GHz range by the USAF Radio 
Solar Telescope Networlj^ (RSTN), the Bumishus patrol telescopes (Institute of 
Applied Physics of Bern University), BEMRAK, the KOSMA telescope, and 
the SST. Above 200 GHz, the radio event is similar to the 2001 April 12 radio 
burst. Indeed: (i) it exhibits a short impulsive phase followed by a time-extended 
component; (ii) the impulsive radio emission is gyrosynchrotron radiation from 
relativistic electrons; and (iii) the gradual phase is optically-thin (8-230 GHz) 
and optically-thick (345 GHz) thermal bremsstrahlung. The main goal of this 
study is to combine the radio observations with soft X-ray, UV/EUV and optical 
data in order to discuss the origins of the thermal radio emission during the 
gradual phase. The paper is organized as follows. Section [2] describes the radio 
observations and flux calibration of the SST, KOSMA and BEMRAK measure- 
ments. The results are discussed in Section [3} and conclusions are drawn in the 
final Section. 



2. Observations and Data Analysis 

The 2003 October 27 radio event at -12:30 UT was associated with a GOES 
M6.7-class soft X-ray (SXR) flare and an Ha sub-flare which occured in active 
region AR 10486 at S17 E25. The flare was well-observed in the UV and EUV, 
from ~12:28 to ~13:00 UT, by the Transition Region and Coronal Explore'^ 



(TRACE; Handy et al, 1999 ). Unfortunately, only the late part of the flare (after 
^12:46 UT) was observed in the hard X-ray domain (>20 keV), by the Reuve n 
Ramaty High Energy Solar Spectroscopic Image^ (RHESSI; Lin et al., 2002). 



Radio observations were available over a wide spectral domain from 10 kHz up to 
405 GHz. Figure [l] shows the time evolution of the GOES 1-8 A SXR flux and of 
the radio flux densitief]^ at selected frequencies in the 1-400 GHz domain. Total 
flux densities have been obtained by RSTN at eight frequencies in the 0.245- 
15 GHz range and by the Bumishus patrol telescopes at 8.4, 11.8, 19.6, 35, and 
50 GHz (until 12:57 UT). Measurements at 212 and 405 GHz were made by the 
SST. KOSMA and BEMRAK provided measurements at 230 and 345 GHz, and 
at 210 GHz, respectively, until 13:00 UT. 

The radio event exhibits two phases: an impulsive burst which lasts for about 
3 min (^12:29-12:32 UT) and a slowly varying and time-extended emission 



^http 
http 



/ /www. ngdc.noaa.gov/stp/solar/solarradio. html 
/ /trace. Imsa Lcom/ 



http://hesperia.gsfc.nasa.gov/hessil 

sfu = 10-22 W m-2 Hz-i 
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Figure 1. Time histories of X-ray and radio emissions during the 2003 October 27 flare at 
~12:30 UT. From top to bottom: 1-8 A soft X-ray flux (background subtracted) from the 
GOES-12 X-ray Sensor; radio flux densities at 1.415 GHz from RSTN; at 19.6 GHz from 
the Bumishus telescopes (the blue curve shows the time extended phase with an enlarged 
flux density scale); at 210 GHz from BEMRAK (red), 230 GHz from KOSMA (blue), and at 
212 GHz from SST (black); and at 345 GHz from KOSMA. 
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which lasts up to about 13:15 UT at 212 GHz. It should be noted that: (i) al- 
though this time-extended emission is more pronounced and larger than the 
impulsive burst in the 200-400 GHz range, it is also observed at all other frequen- 
cies above 1.415 GHz, except at 50 GHz due to strong atmospheric absorption; 
and (ii) at 212 GHz, the slowly-varying emission starts slightly before the impul- 
sive rise, simultaneously with the SXR emission. Such a time sequence is similar 
to that observed in other events with >200 GHz emission (Liithi, Magun, and 



Miller, |2004[ Trottet et a/., |2008[ Gimenez de Castro et aL, |2009[ ). At microwave 
frequencies, two impulsive bursts, Bl and B2, with maxima at ^-^12:35:49 UT and 
^12:42:20 UT, respectively, are superimposed on the slowly-varying component. 

Figure [T] shows that no significant emission is measured at 1.415 GHz. Al- 
though radio activity was observed in the decimeter-hectometer range (see Radio 
Monitoring, these radio emissions are not associated with AR 10486, so the 
flare under study is radio-silent at frequencies below ^2 GHz. Indeed: 

- There is a permanent noise storm which exhibits an enhancement around 
12:30 UT, i.e., close to the maximum of the impulsive phase of the event 
under study. Radio imaging by the Nangay Radioheliograph (NRH) in the 
150-450 MHz range shows that this noise storm activity is spatially asso- 
ciated with AR 10488 (at N09 E12), not with AR 10486. A decameter- 
hectometer Type III storm is also associated with this noise storm, as is 



sometimes the case (see Elgar0y, 1977). 



- Between 12:30 and 14:00 UT, there are some decimeter-hectometer Type III 
bursts, indicating that some electron beams accelerated in the low corona 
reach interplanetary space. Here, again, NRH images show that these 
Type His are not associated with the flare under study, but rather with 
activity from AR 10484 (at N07 W46). 

The 2003 October 27 flare thus belongs to the class of "confined microwave 



events" (Klein et al., 2011 ), where microwave-emitting electrons remain confined 



in the low corona. This is consistent with the analysis of the magnetic field 



topology of this flare by Luoni et al. { 2007 ) , which indicates that energy release 



occurred in compact loops at a magnetic null very low in the corona (^3 Mm 
above the photosphere). 

2.1. Flux Densities at 212 and 405 GHz 

The SST observed the 2003 October 27 event with flve independent total-power 
receivers: channels 1-4 at 212 GHz and channel 5 at 405 GHz. Figure [2] shows 
the 50%-level contours for beams 2-5 while SST was tracking AR 10486, overlaid 
on a TRACE 1600 A image obtained at ~12:46 UT, close to the maximum of 
the time-extended phase. Beam 1 is outside the field of view of the Figure. At 



,http://secchirh. obspm.fr/survey. php?dayofyear=20031027&composite=l and 
|http://secchirh.obspm.fr/survey.php7dayotyear=20031027&^composite=2| 
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Figure 2. Configuration of the SST beams during the fiare. The circles and the ellipse rep- 
resent the 50% contour (HPBW) of beams 2-5, overlaid on a TRACE 1600 A image taken 
at 12:46:21 UT. The numbers indicate the beam. The green triangle represents the region 
where the 212 and 405 GHz flare sources may be located without being detected by beam 2 
at 212 GHz or by beam 5 at 405 GHz (see text). 

212 GHz, beams 1-4 may be approximated by circular Gaussians with a half- 
power beam- width (HPBW) of ~4'. The 405 GHz channel was, for this period, 
oblate, and is better represented by an elliptical Gaussian with a HPBW of 
~2' X 4'. 

Due to the low solar elevation angle 9ei ~ 35° at the time of the flare, along 
with the rather large zenith opacity, the antenna temperature time profiles at 
212 and 405 GHz are quite similar, dominated by atmospheric fluctuations as 
illustrated in Figure |3] However, we note that the signal from beam 3 shows 
additional temporal features compared to those in beams 2, 4 and 5. In partic- 
ular, we observe a rapid variation at 12:30 UT followed by a gradual increase 
between 12:32 and 13:00 UT, which correspond, respectively, to the impulsive 
and extended phases of the flare, shown in Figure [ij Therefore, the flare was 
detected mainly within beam 3, with negligible contributions in beams 2, 4 and 
5. This is consistent with the beam positions relative to the UV flaring structures, 
as shown in Figure [2j 

Since the signal from beam 4 is very noisy, including spiky interferences, we 
instead consider that the signal from beam 2 is the one which is characteristic 
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Figure 3. Antenna temperatures observed by the different SST beams at 212 GHz (beam 2: 
red; beam 3: thick black; beam 4: dashed black) and at 405 GHz (beam 5: blue). Arrows 
indicate the maxima of the flare impulsive and extended phases at 12:30 UT and 12:45 UT, 
respectively. 



of the emission of the whole active region. Figure [3] shows that, before the flare, 
the signals from beams 2 and 3 are identical. This indicates that the pre-flare 
(background) signal does not depend critically on the position of the whole active 
region within the beam. The signal of beam 2 has thus been subtracted from 
that of beam 3 in order to suppress variations due to atmospheric fluctuations 
and to obtain the excess antenna temperature due to the flare. 
The flux density at 212 GHz, S'212, is then given by: 

<S'2i2 = ^v^exp(rz/ sin (6*01)) , (1) 
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where /cb is the Boltzmann constant, Aa is the antenna surface area, r]a the 
aperture efRciency, and the zenith opacity of the Earth's atmosphere. We 
estimated T2(212) ~ 0-46 and T2(4g5) « 1.9 at 212 and 405 GHz, respectively, 
using the method described in Melo et al. (2005). For the present observation, 
the antenna aperture efficiency rja is 20% at 212 GHz. The main sources of 
uncertainty are temporal fluctuations of and the aperture efficiency, which 
result in a total uncertainty of 20%. Figure [T] shows the 212 GHz flux density as 
computed using Equation Q. 

In Equation ([T]) , we have assumed that the source is located over the center of 
beam 3. In order to account for the actual source location, we need to multiply 
S212 by the factor exp{9lff/2a^). Here, ^off is the angular distance between the 
center positions of the source and of beam 3, and cr represents the standard 
deviation of the convolution of a circular Gaussian source by the beam. We 
estimated that the maximum flux density at 405 GHz should have been ^125 sfu 
by extrapolating the increasing 212-345 GHz flux spectrum up to 405 GHz 
(see Figure [I]). However, the event was detected only in beam 3 at 212 GHz, 
and was not detected in beam 5 at 405 GHz within the r.m.s. noise of 30- 
40 sfu. We use these two facts to constrain the source size and location, and to 
estimate 6'ofi. To do so, we further supposed that the center of the 405 GHz- 
emitting source should be near the maximum of the SXR emission observed by 
RHESSI (see Section [2?3]) and between the UV ribbons observed by TRACE (see 
Figure [2]). We then computed the emission that should be detected at 405 GHz 
by convolving a circular Gaussian source of 125 sfu with beam 5. We found that 
sources with sizes ranging from 10" to 70", with centers located within the green 
triangle in Figure[2] fulfill the above two requirements, i.e., no burst detection in 
beam 5 at 405 GHz or in beam 2 at 212 GHz. The correction factor for 5*212 in 
Equation ([I]) is thus between 25% and 40%, including a 20" indetermination in 
antenna pointing. At 212 GHz, the maximum uncertainties on the flux density 
are thus taken as -30% and +40%. 



2.2. Flux Densities at 210, 230 and 345 GHz 

On 2003 October 27, the KOSMA target was AR 10486, where the flare un- 
der study was located. KOSMA uses a tillable subreflector which periodically 
deflects the antenna beams by 6' in elevation, from the flaring active region 
to a quiet reference region on the Sun. The "1 second on / 0.1 second off" 
observation cycle allows for compensation of atmospheric attenuation changes 
with time scales greater than 1 s (Liithi, Magun, and Miller, 2004). The absolute 
flux density calibration was obtained from hourly observations of the sky. Unfor- 
tunately, before 12:31 UT the subreflector was wobbling the wrong way, staying 
10 s on the quiet reference region and 0.1 s on AR 10486. With such a short 
sampling of the active region, we could determine neither the excess flux of the 
pre-flare active region with respect to the quiet reference flux, nor the onset of 
the flare itself. Signiflcant flux density excesses at 210, 230 and 345 GHz relative 
to the quiet reference region were thus measured only after 12:29:56 UT, i.e., 
slightly before the maximum of the impulsive phase at 212 GHz (see Figure [T]). 

In order to estimate the net flux density excess due to the flare alone, we 
proceeded as follows. We flrst assumed that the flare emission started at the 
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same time at 210, 230 and 345 GHz as at 212 GHz. This is generally the case for 



events detected in the 200-400 GHz range (e.g., Kaufmann et al., 2004 Liithi 



Magun, and Miller, 2004 Rauhn et al, 2004 Kaufmann et al., 2009). The growth 



of the time-extended emission at 212 GHz is well-fit by a straight line between 
12:32 and 12:35 UT. By extrapolating this straight line back to the time to where 
it crosses the pre-event flux density level, we get an estimate of the starting time 
of the burst: to = 12:38:37 UT ± 5 s. At 210 and 230 GHz, straight hues have 
been fit to the same time interval as at 212 GHz and extrapolated back to tg, 
providing the actual pre-flare flux densities at 210 and 230 GHz. These are found 
to be 36 ±6 and 47 ±7 sfu at 210 and 230 GHz, respectively. These uncertainties 
are estimated by taking into acount the uncertainty on to, the uncertainties on 
the fit parameters, and the amplitude of the fluctuations of the measured flux 
density. Figure [l] shows that the derived flare excesses at 210 and 230 GHz agree 
well with that measured at 212 GHz, justifying the procedure described above. 

At 345 GHz, the early increase of the time-extended emission cannot be fit by 
a single straight line. We obtained upper and lower limits for the contribution 
of the active region before the flare by fitting straight lines to the 345 GHz flux 
density excess measured between 12:32 and 12:32:34 UT, and between 12:32:34 
and 12:33:50 UT, respectively, and extrapolating back to to- The contribution 
of the active region is thus found to be between 74 ± 8 and 98 ± 10 sfu. For the 
following Sections, we adopt the mean value, 86 ± 14 sfu. 

Following Liithi, Magun, and Miller ( |2004[ ), we adopt maximum uncertainties 
of -30% and +A0% at 210 and 230 GHz. An absolute error of ±14 sfu is added 
at 345 GHz (see previous paragraph). Figure [T] shows that, during the time- 
extended phase, the time evolution of the 210-230 GHz flux and of the 345 GHz 
fluxes are different, with the flux density at 345 GHz being systematically larger 
(about 2x after 12:40 UT) than that at 210-230 GHz, even when using the 
upper limit of the pre-flare active region contribution. As the above uncertainties 
constitute conservative upper limits of the flux accuracy, we consider that the 
difference between the 210-230 GHz and 345 GHz flux densities is significant 
despite being only slightly larger than the uncertainties. 

2.3. RHESSI spectral and imaging observations 

We used RHESSI X-ray spectra and images to quantitatively characterize the hot 
(>10 MK) thermal plasma in the 2003 October 27 flare. While the GOES SXR 
data provide some information about the flare plasma temperature and emission 



measure (see Section 3.2.1 ), the broadband response over only two channels (1- 
8 and 0.5-4 A, corresponding to ^1.6-12 and ^3.1-25 keV, respectively) thus 
yields only a single, average measurement with no spatial information. On the 
other hand, in the X-ray range (~3-100 keV), RHESSI's -1 keV FWHM spectral 
resolution enables precise measurements of the hottest plasma, which can then 
be combined with the GOES observations to also measure cooler components 
and to thus more accurately determine the temperature distribution in the flare. 
RHESSI also provides spatial information for the hot plasma, with an angular 
resolution as good as ~2". 

RHESSI observations of the 2003 October 27 flare are available only in the 
late decay, from ^12:46 to ^-^13:04 UT. Because of the limited observing period 
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and the rapidly-decaying hard X-ray flux (the incident ^6-25 keV flux decreases 
by about two orders of magnitude over these ~18 minutes), we analyzed only a 
single time interval, from 12:47 to 12:49 UT. The 3-100 keV spatially-integrated 
spectrum was accumulated over this 2-minute period using the single best de- 
tector, G4, with 1/3-keV energy bins (the instrument channel width; Smith 
et al., 2002 ), and the non-solar background was then subtracted (cf. Caspi, [20T0 ); 
above ~3 keV, the solar pre-flare background is negligible compared to the flare 
emission. The spectrum decreases steeply with energy (e- folding of ^2 keV), 
with no appreciable flare emission above ~33 keV in the analyzed period. 

For precise spectroscopy, we employed the low-energy instrument response 
calibration improvements described by Caspi (20101. Following the method of 
Caspi and Lin ( |2010 ), the ~4. 67-33 keV spectrum was forward-flt with a photon 
model consisting of two isothermal continua (a single isothermal did not yield 
an acceptable fit) and two Gaussian features (representing unresolved excitation 
lines of highly- ionized Fe and Ni). The spectrum was well-fit by this model 
(reduced ~ 0.71) and yielded, for each of the two thermal components, a best- 
fit temperature and emission measure of Ti w 24 MK, EMi « 0.24 x 10"'^ cm^'^, 
Ta « 16 MK, and EM2 ~ 1.7 x 10''^ cm^^ respectively. 

An image at 6-30 keV was accumulated for the same time interval, using the 
CLEAN image reconstruction algorithm with uniform weighting (Hurford et al, 



2002) and grids 3, 4, 5, 6, 8, and 9 (grid 7 was excluded due to its ^20 keV 



low-energy threshold) to achieve a ~10" spatial resolution (see Section 3.2.1 ). Al- 
though it is not possible to image the two thermal components separately, we can 
nevertheless derive their centroid positions relative to the combined emission (cf. 
Caspi and Lin, 20101. From the spatially-integrated spectral model, the ^-^24 MK 
source was found to contribute ~45%, ~63%, and ~90% of the total count fiux at 
6.3-7.3, 9-12, and 17-18 keV, respectively, and imaging shows that the centroid 
position of the emission in each energy band is correlated with this fractional 
contribution. From this, we extrapolated the centroid positions of the ~24 and 
~16 MK thermal plasmas and found them to be separated by ~11" ± 3.2". The 
6-30 keV source is dominated (^S0%) by counts from the ~24 MK component, 
and encompasses a projected area (cf. Caspi, 2010) of ~280-310 square arcsec 
within the 50% contour, and ~560-660 square arcsec within the 30% contour 
(both values have been corrected for broadening by the instrument point-spread 
function). Although the ^16 MK component contributes only weakly to the 
imaged emission, the morphology and centroid separation suggests that it is of 
a similar size. 



3. Results and Discussion 
3.1. Impulsive Phase 

Figure |4] displays the radio spectrum observed in the 1-400 GHz range around 
the maximum of the impulsive phase at 212 GHz. This spectrum is reminiscent 
of the usual gyrosynchrotron emission from non-thermal electrons, commonly 
observed in the microwave domain. It peaks around 10 GHz, has a sharp low 
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Figure 4. The radio spectrum near the maximum of the 212 GHz emission during the im- 
pulsive phase of the 2003 October 27 event. The solid line is a power-law fit to the decreasing 
part of the radio spectrum, with spectral index a = —1.2. 



frequency cutoff between 1.4f 5 and 2.695 GHz, and has a positive spectral index 
of ~4 below 5 GHz. This is steeper than the index of 2.5-2.9 predicted by 
gyrosynchrotron self-absorption. Such steep spectra are generally ascribed to the 
suppression of microwave emission by the ambient medium (the Razin effect). In 
magnetic fields of a few hundred Gauss, the Razin effect becomes significant at a 
few GHz for ambient densities of at least a few lO^*' cm^'^. For stronger magnetic 
fields, higher ambient densities are required for the same level of suppression. 
Such high densities inferred for the microwave-emitting region further support 
the statement that the 2003 October 27 flare is a "confined microwave event" 
(see Section [2]). 

The optically-thin, decreasing part of the spectrum above 15 GHz is well- 
represented by a power law, as shown by the solid line in Figure[4] The power-law 
spectral index a is about —1.2 around the maximum of the impulsive phase at 
212 GHz. This optically-thin radiation is primarily emitted by MeV electrons. 
Considering the ultra-relativistic case as a gross approximation, the spectral 
index 6 of the instantaneous distribution of radio-emitting electrons is 5 = 2q; — 1 
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(Ginzburg and Syrovatskii, 1965[ p thus 6 ~ —3.4 at the maximum of the impul- 



sive phase. Such a hard electron spectrum is characteristic of gyrosynchrotron 
events observed up to millimeter-submillimeter wavelengths (e.g., Trottet et al., 
[2002{ Liithi , Magun, and Miller, ,2004, Raulin et al, |2004( Gimenez de Castro 
et al., 2009 and references therein). 

Figure |4 indicates that, within the large uncertainties, the emission at 210- 
230 and 345 GHz appears to be the high-frequency part of the gyrosynchrotron 
spectrum observed in the microwave domain. Such a continuation of the mi- 
crowave spectrum to the 300-400 GHz domain has been previously observed, 
although for larger, GOES X-class flares. 

3.2. Gradual Phase 

As stated in Section [2] the slowly- varying and time-extended ( "gradual" ) radio 
emission starts to rise close to the onset of the SXR emission measured by GOES 
(see Figure [T]). Superimposed on this slowly-varying component, two impulsive 
bursts, labeled Bl and B2 on Figure [l] are observed at frequencies up to 35 GHz 
but without obvious counterparts in the 200-400 GHz domain. Figure [s] shows 
the mean radio spectrum (asterisks) observed during Bl and B2 and during two 
later time intervals. For each of the four time intervals shown, the dashed curve 
indicates the mean bremsstrahlung spectrum expected from a single isothermal 
source of assumed 40" diameter with temperature T„ and emission measure EM„ 



as derived from the GOES SXR fluxes (see Section 3.2.1). The main points to 
be drawn from Figure [T] and Figure |5] are: 

- During the two last intervals, when no impulsive bursts are present, the 
observed radio spectrum is almost flat between ~8 and 230 GHz, with flux 
densities ^^3 to 4 times greater than those expected from the single thermal 
source derived from GOES SXRs. 

- During bursts Bl and B2, a gyrosynchrotron component is superimposed 
on a thermal bremsstrahlung spectrum. Bl has no counterpart in the mm- 
submm range, but the >200 GHz radiation exhibits an additional increase 
during B2. 

- Throughout the gradual phase, the 345 GHz flux density is significantly 



higher than that measured in the 210-230 GHz range (see Section 2.2). 
We discuss each of these points in the following Sections. 
3.2.1. The Flat Radio Spectrum 

The flat spectrum observed between ^8 and 230 GHz throughout the gradual 
phase is reminiscent of optically-thin thermal bremsstrahlung emission. Below 



^Ginzburg and Syrovatskii 1 1965 \ define «, 5 > 0; we have modified their equation to eliminate 
this assumption, taking a — > —a and 5 — > —5. 
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Figure 5. The mean radio spectrum observed at different times during the time-extended 
(gradual) phase of the 2003 October 27 flare (asterisks with errors bars). From top to bottom, 
the spectra correspond to bursts Bl and B2 (see Figure[ljl, and to the plateau and decay of the 
210-230 GHz emission, respectively. The dashed lines indicate the radio spectrum expected 
purely from an isothermal plasma with temperature and emission measure derived from the 
GOES SXR measurements, assuming a circular source of 40" diameter, while the solid lines 
indicate the spectrum expected from a two-temperature model (see text). 



~8 GHz, the flux density decreases with decreasing frequency, indicating that 
the radio source graduafly becomes opticafly-thick at low frequencies. To char- 
acterize the thermal plasma seen by GOES, its temperature Tg and emission 
measure EMg were derived from the GOES SXR fluxes, assuming an isothermal 
source, using the CHIANTI database (Dere et al, 1997 Young, Landi, and 
Thomas, 1998) and coronal elemental abundances. Tg peaks at ^-^20 MK just 
after the maximum of the impulsive phase, then decreases to ^12 MK. EM„ 
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initially increases, passes through a maximum of ^3 x 10*^ cm^'^ approximately 
at the end of B2, and then falls to ^lO"*^ cm"'^. 

We computed the radio spectrum expected from the GOES thermal source as 
follows. At time t and frequency v, the brightness temperature Tbr(i, i') of the 
radio emission from an isothermal, homogeneous source at temperature T{t), 
emission measure EM{t), and projected surface area A (normal to the line-of- 
sight) is given by: 

Tb.(t,i^)=T(t)(l-e-^^'^)) , (2) 
where, for T > 0.2 MK, the optical thickness r(i, v) is given by (Dulk, [19851 ): 
K{t, v) EM{t) 



^^^'"^ v'^A T{t)y^' (3) 
K{t, J/) = 9.78 X 10"^ X (24.5 + In T{t) - In i^) (c.g.s. units) . 

The expected radio flux density S'cxp(^, i^) from the source is then given by: 

where c is the speed of light and R is the Sun-Earth distance. We note that for 
optically-thin (t ^ 1) regimes, Equations ^ and Q reduce to: 



2kBK{t, v) EM{t) 



whereby the expected radio flux spectrum is independent of the source area A 
and only weakly dependent on frequency v. 

Using Tg(t) and EM^{t) derived for the GOES thermal plasma, and assuming 
a projected area equivalent to a circle of 40" diameter (the mean acceptable 
source area, per Section |2.1[ see also the discussion below) , we calculated the 
radio flux density Scxp{t, u) expected for this source, shown by the dashed black 
curves in Figure [5j Like Liithi, Magun, and Miller ( |2004[ ), we find a qualitative 



similarity between the time evolutions of the observed and predicted flux densi- 
ties in the 200-400 GHz range (where the emission is optically thin and where no 
impulsive bursts are detected), but the observed values are ~3-4 times higher at 
210-230 GHz and ~7 times higher at 345 GHz. Even using photospheric rather 
than coronal abundances to derive Tg and EMg, the computed radio fluxes are 
still ~2 and ^4 times lower than the observed values at 8-230 and 345 GHz, 
respectively. 



In contrast, Pohjolainen, Valtaoja, and Urpo (1996) reported gradual radio 
bursts at 37 GHz for which the observed and computed flux excesses were com- 
parable. A similar agreement was obtained from case studies at 8.8 and 15.4 GHz 



by, e.g., Kundu et al. (1994), up to 15.4 GHz by Trottet et al. (2002), and up to 



86 GHz by Raulin et al. ( 1999 ). However, we note that in these earlier studies, Tg 



and EMg were derived using the Thomas, Crannell, and Starr ( 1985 ) polynomial 
approximation, which closely resembles the GOES response to photospheric- 
abundance CHIANTI models for temperatures up to ^15 MK (White, Thomas, 
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and Schwartz, 2005). The GOES plasma Hkely resides in coronal loops, not chro- 
mospheric footpoints, and as shown above, the use of photospheric abundances 
results in substantially larger expected radio flux densities than those obtained 
using the more appropriate coronal abundances. Thus, this suggests that, for 
most of the gradual bursts reported in the cited works, the measured flux densi- 
ties in the centimeter-millimeter range are actually higher than those expected 
from the SXR-emitting plasma observed by GOES, similar to our results. 



White and Kundu ( 1992 ) emphasized that the above discrepancy between ob- 
served radio emission and that computed from SXR observations mainly reflects 
that: (i) the assumption of an isothermal coronal plasma is not correct during 
a flare, and (ii) radio observations are more sensitive to cooler material than 
is GOES because, in contrast to SXR intensity, the intensity of optically-thin 
radio emission decreases with increasing temperature. They concluded that a 
0.1-1 MK source would produce a stronger response in the centimeter-millimeter 



domain than would (hotter) plasma observed by GOES. Chertok et al. (1995) 
reached similar conclusions in a joint analysis of GOES observations and radio 
bursts measured at 3-80 GHz. 

Following these authors, we consider that, in addition to the coronal SXR 
source SG observed by GOES, there is an additional, cooler, radio-emitting 
coronal source SR at plasma temperature < Tg with emission measure EM^. 
For simplicity, we assume that SR has the same surface area A, normal to the 
line of sight, as SG and that it is located directly below it. Then, by modifying 
Equation ^ for this geometry, we obtain the brightness temperature of the 
combined source: 

Tta(i, ly) = Tg(i) (l - e-^«(*'^)) + T,{t) (e-^^^*'"^) - e-(^«(*>'')+^'(*>''))) , (6) 

which, together with Equations (|3| and Q, yield the expected coronal flux 
density Sl^^{t, u) for the combined source. If both SG and SR are optically thin 
(rg ^ 1 and Tj ^ 1), Equations ^ and Q reduce to: 

Tb,.(t, u) « Tg(t)Tg(t, v) + T,{t)T,it, v) , 

2fcB fK^{t,u)EMS) , K,{t.v)EAm\ (7) 



:2i?2 \^ Tg(t)l/2 T,(t)l/2 

Tg and EM^ are fixed by the GOES observations; the free parameters Tr, EM^, 
and A can then be determined by fitting the model ^^"(t, v) (Equation Q) to 
the whole range of observed flux densities S'obs(^, v). 

As a first approach, we varied the free parameters of the cool plasma SR 
until a reasonable fit to the data was achieved by eye. We constrained the fit 
parameters by requiring that the cool plasma does not contribute significantly 
to the GOES SXR flux in either energy channel, since that would conflict with 
the previously-derived Tg and £Mg which were held fixed in the model. In this 
way, we determined the best-fit parameters of the cool plasma for 1-minute time 
intervals from 12:33 to 12:57 UT; during Bl and B2, when synchrotron radiation 
adds to the thermal bremsstrahlung up to at least 35 GHz, Tr and EM^ were 
adjusted to account for the observed 210-230 GHz flux densities. 
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Figure [5] illustrates the reasonable agreement of the two-temperature coronal 
source model (solid black line) with the observations (asterisks) for A = 6.7 x 
10^® cm^, which corresponds to a diameter of 40" for a uniform circular source. 
Tr and £Mr, respectively, begin at 1 MK and 7 x 10**^ cm~^ in the first time 
interval (12:33-12:34 UT), reach a maximum of 2.6 MK and 3.5 x 10*^ cm'^ in 
the 12:43-12:45 UT time period (decay of burst B2), and decrease to 1.3 MK and 
1.2 X 10"*^ cm-3 at the end of the analyzed period (12:57 UT). This 40", 1-3 MK 
plasma, with EM^ « EMg, contributes <1% of the GOES SXR flux in either 
channel, as we required. Larger sources, up to 60", yield similar values of 
and EMj., but smaller sources yield higher temperatures and emission measures, 
resulting in a non-negligible predicted contribution of the SR source to the GOES 
flux in at least the low-energy channel. A 30" source, for example, with best-fit 
T, « 4.6 MK and EM, « 5 x 10^^ cm'^ « 2 x EMg, would contribute only 
-0.5% of the SXR flux at 0.5-4 A but a significant -16% at 1-8 A, which 
violates our initial constraints. Thus, we adopt a coronal radio source diameter 
D, = 40" ± 10". 

While this simple two-temperature model prediction of the radio spectrum 
yields reasonable agreement with the observations, the RHESSI data, when 
available, provide distinct evidence that the true temperature distribution of 
the thermal flare plasma is more complex. For the specific time interval of 



12:47-12:49 UT (see Section [2^, the RHESSI -3-33 keV spectrum is well 
characterized by two hot, isothermal components with Ti « 24 MK, EMi « 
0.24x10^^ cm-^ T2 w 16 MK, and £^2 - 1.7x10^^^ cm^^. Because of the broad- 
band GOES response, the —24 MK source will contribute non- negligibly to the 
measured SXR flux, which will therefore affect the radio spectrum predicted from 
the GOES measurements. In particular, the two thermal components observed 
by RHESSI contribute a combined 85% and 98% of the GOES SXR flux ob- 
served in the low- and high-energy channels, respectively. The remaining GOES 
fluxes are consistent with an additional "warm" plasma, with T3 « 6.7 MK and 
EM^ « 1.0 X 10^^ cm~^. However, the red dashed curve in Figure [s] shows that 
the radio emission predicted from this three-temperature (Ti, T2, T3) plasma is 
still below the observed spectrum. Thus, a fourth, cool plasma, with T4 « 2.0 MK 
and EM4 « 3.0 x lO''^ cm~^, is still needed to account for the observed radio 
flux densities (with a diameter D, = 40" ± 10", as before); this cool plasma 
contributes negligibly (<1%) to the GOES fluxes, as required. 

Figure [5] shows that the four-temperature model prediction (solid red line in 
the 12:47-12:49 UT panel) is nearly identical to that of our previous, simple, 
two-temperature model (solid black line). Thus, in the absence of RHESSI data, 
a simple two-temperature scenario nevertheless constitutes a viable working 
model, representative of the radio spectrum from a multi-temperature coronal 
plasma. However, we note that the derived T and EM, and hence the pre- 
dicted radio contribution, of the cool source do vary between the two models, so 
additional measurements within the temperature distribution are desirable for 
optimal accuracy, when possible. 

TRACE observations at 195 and 1600 A were made between 12:28 and 
13:00 UT, with a few observations at 284 A also made after 12:41 UT. Images 
at 195 A (and at 284 A) show the geometry of the flaring coronal plasma, while 
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t igure 6. TRACE images at 195 A ( reverse grayscale) and 1600 A (green contours) at four 
instants during the 2003 October 27 flare; the TRACE alignment was corrected by [+6", 
-7.5"]. The 12:30:45 UT image (top left) occurs during the impulsive phase; r4, r 5. r7 a nd r8 
refer to the 1600 A flare ribbons/kernels as identified in Figure 4 of Luoni et al. (20071, and 
the gray b oxes indicate the areas over which the kernel intensities have been computed (see 
Sec. |3.2."3| . The three subsequent images, which exhibit two loop systems marked LI (south) 
and L2 (north), occur during the gradual phase. RHESSI 6-30 keV image contours (30%, 
50%, 70%, 90%) during 12:47-12:49 UT are overlaid on the final image (bottom right), along 
with the derived centroid positions of the two thermal plasmas inferred from the simultaneous, 
spatially-integrated RHESSI spectrum (see Sec. |2.3| l. For comparison, the dashed black circle 
of 40" diameter represents the model coronal radio source, arbitrarily centered on the centroid 
of the ~16 MK plasma. 



the 1600 A images display the chromospheric footpoints/ribbons, as illustrated 
in Figure [g] at four instants during the event. The bright 1600 A kernels are 
marked r4, r5, r7, and r8, as in Figure 4 of Luoni et al. (2007); similar kernels 
are also seen in Ha (see Figure 3, right panel, of Luoni et al., 2007). The image 
at 12:30:45 UT, characteristic of the impulsive phase, shows only faint 195 A 
emission directly atop the 1600 A kernels. During the gradual phase, Figure [g] 
shows two bright 195 A loop systems, marked LI and L2, which connect the 
1600 A kernels. The 195 A intensity, integrated over the area covered by these 
loops, increases from the beginning of the gradual phase up to at least 12:40 UT; 
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after that time, the images are contaminated by numerous spikes that prevented 
a rehable measure of the Hght curve. 

The 195 A passband is sensitive to spectral hue emission from both Fe XII 
(peak contribution at ~ 1.1-1. 7 MK, cf. Mazzotta et al, |1998| ) and Fe xxiv 
(~13-35 MK), and the 195 A images wih thus contain some contribution from 
the hot ~24 and ~16 MK plasmas seen by RHESSI. However, per unit emission 
measure, the se nsitivit y is ^10 times stronger at the low (Fe xii) temperatures 



(cf. Handy et al., 1999 1, and thus the total 195 A emission in LI and L2 should be 
vastly dominated by cool plasma. This is confirmed by the similar loop structures 
observed at 284 A (Fe xv, ~1.8-2.6 MK), and suggests that LI and L2 are indeed 
largely representative of the derived 1-3 MK source considered in our working 
model, and therefore of the coronal thermal radio-emitting region. The increase 
of the 195 A intensity up to at least 12:40 UT is thus consistent with our model, 
which predicts an increase of Tj. and EM,- until 12:43-12:45 UT. 

The spatial distribution of the RHESSI-observed 6-30 keV emission during 
12:47-12:49 UT, where we applied the four-temperature thermal model, is also 
shown in Figure [6] (bottom right panel) . The X-ray source has a similar shape 
as and projects onto the LI loop system. The projected centroid positions of 



the ^24 and ~16 MK plasmas (see Section 2.3) are located near the top of the 



brightest loops within LI, as would be expected for the hottest plasma (cf. Caspi 



and Lin, 2010). The projected (and deconvolved) surface area of the RHESSI 
X-ray source, ~300 and ~600 square arcsec within the 50% and 30% contours, 
respectively, is significantly smaller than the ^1250 square arcsec area of the 
radio source obtained from the model. Since the X-ray source produces only 
^20% of the radio flux (see Figure [5|, its size is not likely representative of the 
primary radio source, suggesting that the 4-T model assumption of identical sizes 
for the hot (>16 MK), warm (~6-7 MK), and cool (~l-3 MK) sources may not 
be valid. However, this does not affect the results in the optically-thin regime, 
which are independent of source size (cf. Equation ([t])). Although the 195 A- 
emitting region (bounded by a ^60" x 50" rectangle) is very inhomogenous, with 
a diffuse edge that is difficult to define precisely, we estimate that the projected 
surface area is between ~800 and ~1150 square arcsec. While the areas of the 
X-ray and EUV sources are estimated using dissimilar methods, the difference in 
their areas is significant, and it is noteworthy that the X-ray source has negligible 
emission overlying the L2 system, which can largely account for its smaller size. 
The 195 A source area is, on the other hand, in reasonable agreement with the 
^1250 square arcsec area of the model radio source, consistent with the model 
finding that the radio emission is primarily from the ~l-3 MK plasma observed 
by TRACE. Unfortunately, the present data do not provide enough constraints to 
consider a more sophisticated model which accounts for inhomogeneities within 
the EUV- and X-ray-emitting regions; observations from UV/EUV and X-ray 
imagers and spectrometers onboard the Solar Dynamics Observatory and Hinode 
satellites may provide the necessary information to develop such complex and 
realistic models in the future. 

In summary, the 1-230 GHz radio spectrum observed during the gradual 
phase is consistent with thermal bremsstrahlung radiated by a coronal source 
of ^-^40" equivalent diameter, filled with hot, ^7-24 MK plasma observed by 
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GOES and RHESSI and with cooler, -1-3 MK plasma observed by TRACE. 
The similarity between the time evolution of the SXR and radio emissions, par- 
ticularly in the 200-400 GHz range (see Figure [T]), implies that these emissions 
likely share a common origin. The increase of EMg and EM^ from the beginning 
to the maximum of the gradual phase suggests that this common origin could 
be chromospheric evaporation. 

Since the radio corona is optically-thin above —8 GHz, the observed flux 
density S'obs(i') at frequency > 8 GHz can be written as the sum of the emission 
from the corona, S^^Ki^), and from the chromosphere, ^'^{^^(z/). Similarly, the flux 
density 5*0 before the flare is thus S*™'' + Sq^'-' . Then, the excess flux density due 
to the flare, S'f, is given by: 

5f(^) = 5oTs('^) - Sri'^) + S^ti'^) - . (8) 

The observation of optically-thin bremstrahlung during the gradual phase leads 
to the following statements: (i) the exce ss flux seen at 345 GHz and the enhanced 



1600 A and Ha emissions (see Section |3.2.3[ below) provide evidence that the 



chromosphere responded markedly to flare energy deposition, but (ii) Equa- 
tion ([s]) implies that, since the observed flux densities are well-modeled entirely 
by coronal sources (per above), 5'obs('^) ~ Sq^'^{i^)] that is, the chromosphere did 
not contribute substantially to the time-extended radio emission over a large 
frequency range Au « 10-230 GHz. This range may vary from flare to flare, and 
even during a single flare, such as was the case during the 2001 April 12 event 



(Liithi, Magun, and Miller, |2004| ) for which Av varies from ~10-230 GHz at the 
maximum of the thermal phase to ~10-100 GHz in the late decay. Variations 
of Av during the gradual phase likely reflect the dynamics of the chromospheric 
response to flare energy deposition. 

3.2.2. The Bursts Bl and B2 

Although a detailed study of bursts Bl and B2 is not of interest here, it is 
necessary to identify the origin and the nature of these bursts to understand the 
relationship between the radio, UV and Ha emission discussed in Section f3. 2. 3[ 
below. 

Bl and the smaller burst just prior to it (see Figure [T]) are both observed 
between 5 and 35 GHz in time-coincidence with Type HI electron beams recorded 
by NRH. NRH imaging in the metric domain shows that these Type Ills arise 
from AR 10484 at N07 W46. Furthermore, high cadence Ha images, obtained 



with the Ha Solar Telescope for Argentina (HASTA; Bagala et ai, 1999 Fernan- 



dez Borda et al., 2002), also show simultaneous brightenings in the same active 
region. We thus conclude that Bl is not part of radio emission arising from the 
2003 October 27 flare under study. 

In contrast, B2, which is also observed in the 5-35 GHz range, does arise 
from AR 10486. Indeed, Figure [7] shows that B2 is coincident with an impulsive 
increase of the four 1600 A kernels marked r4, r5, r7 and r8 in Figure 6 The 
time-derivative of the GOES 1-8 A SXR flux, which may be considered as a 
proxy for non-thermal hard X-ray emission (cf. the "Neupert effect;" Neupert, 
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1968 Dennis and Zarro, 1993), also exhibits a bump coincident with B2 (see 
Figure]?]). As both non-thermal X-rays and the 1600 A UV continuum are widely 
believed to result from direct particle injection into the chromosphere (Fletcher 



and Hudson, 2001 Coyner and Alexander, 2009), B2 is thus the signature of 
a new episode of particle acceleration in the AR 10486 magnetic field. After 
subtraction of the thermal radio spectrum expected from the model described in 



Section 3.2.1 the remaining B2 spectrum is a typical gyrosynchrotron spectrum 
with a maximum at ^10 GHz and a power-law spectral index of approximately 
—2 above that frequency. For such an index, the expected flux density above 
200 GHz is <1 sfu at the maximum of B2. This is consistent with the lack of 
detection of synchrotron emission from B2 at these high frequencies. However, 
the >200 GHz emission shows a significant increase from the beginning of B2 
until slightly after its maximum. This is an indication that the new energy 
release associated with B2 is not only a signature of a new episode of particle 
acceleration, but that it also produced a significant increase of the thermal radio 
emission measured during the gradual phase. 

3.2.3. The flux density increase between 210-230 and 345 GHz 



Luoni et al. (2007) noted that the 195 A-emitting loops shown in Figure |6] with 
footpoint regions r4, r5, r7, and r8, extended above the inversion lines separat- 
ing an elongated region of negative polarity from positive polarities northward 
and southward. Similar kernels, which we label H4, H5, H7, and H8, are also 
seen in the Ha line center. [Luoni et al.\ suggested that the flare energy release 
resulted from magnetic reconnection at a coronal null located at a height of 
3.1 Mm, slightly to the east of r7 (see their Figure 6). In order to investigate the 
chromospheric response to this energy release, we analyzed the time evolution 
of the 1600 A and Ha kernels. For this, we computed the mean intensities per 
pixel over the boxes shown in Figure [6) Figure [7] shows the results for both the 
1600 A kernels (black histograms) and the Ha kernels (red histograms). The 
bottom panel shows the excess flux at 345 GHz with respect to the 210 GHz 
flux ((55345 = 5*345 — 52io)- The examination of Figure [t] suggests the following: 

- Kernels r4 and r5 show similar time behaviors. Indeed, they both exhibit 
a short peak corresponding to the radio impulsive phase, another peak 
corresponding to the radio burst B2, and similar time-decays. H4 and H5 
also evolve similarly, but unlike for r4 and r5, the impulsive phase and 
B2 are not seen in their time proflles. We emphasize that the H4 and H5 
intensities remain roughly constant from ^-^12:42 UT until the end of the 
analyzed period. 

- There is a close similarity between the r7 and r8 light curves. The marked 
difference compared to r4 and r5 is that the impulsive phase is not clearly 
observed, and their decay times are slightly shorter. H7 and H8 evolve 
together, but unlike H4 and H5, they decay after ~12:43 UT. It should be 
noted that the region in which kernel r7 (H7) is located is highly dynamic. 
The elongated shape of this kernel is indicative of the photospheric trace of 
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Figure 7. Prom bottom to top: time evolution of the excess flux density at 345 GHz with 
respect to that of the flux density at 210 GHz; and time evolution of the four 1600 A kernels 
marked r4, r5, r7 and r8 on Figurem](black curves) and of the corresponding Ha kernels H4, H5, 
H7 and H8 (red curves). At 1600^, the pre-flare intensity of the kernel has been subtracted. 
The intensity excess of a given Ha kernel / — /^ef is given with respect to that of a quiet area 
/ref- The blue curves overplotted on the time histories of kernels r7 (H7) and r8 (H8) are, 
respectively, the time evolutions of the 19.6 GHz radio emission and of the time-derivative of 
the GOES 1-8 A SXR flux. 



the fan associated with the null point found in Luoni et al. (2007). Magnetic 
reconnection at the null point would inject energy into the chromosphere 
at the footpoints of loops anchored in r7, but at the same time, would 
induce a disturbance at the section of a long, curved filament lying along 
the inversion line located in its close vicinity (see, e.g.. Figure 2 in Mandrini 
et al., 2006 1. The temporal evolution of the pixels included in the box shown 
in Figure |6 may thus reflect not only the chromospheric response at the loop 
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footpoints but also the nearby filament activity. Nevertheless, we computed 
the r7 intensity for different sub-boxes within the box shown in Figure |6] 
and found that the shapes of the time profiles are not substantially different 
from the one shown in Figure [7] 

The above observational findings therefore support the following statements: 



After reconnection and energy release at the magnetic null point, the two 
sets of 195 A loops LI and L2 have footpoints at r4 (H4) and r5 (H5) and 
at r7 (H7) and r8 (H8), respectively (see Figure [6]). Taking into account 
the magnetic field evolution (e.g., emergence of a new bipole into the pre- 



existing field configuration), Luoni et al. (2007) proposed that magnetic 
reconnection at the null point would proceed in such a way that field lines 
connecting r4 to r7 and r5 to r8 would correspond to LI and L2, contrary 
to what our analysis of the temporal evolution of the flare kernels show. 
Note that it is very difficult to clearly separate these two loops in the 
horizontal "Y-shape" of the TRACE 195 A brightening. However, compu- 
tation of the magnetic field topology is only an indication of the magnetic 
field connectivity and does not tell in which direction the reconnection 
process proceeds; this is inferred by analyzing the temporal evolution of the 
associated phenomena. Furthermore, magnetic reconnection may proceed 
in one direction and, afterwards, in the reverse direction in consecutive 



events, as was shown by Gofi' et al. (20071. Therefore, we conclude that the 



magnetic field topology computed by Luoni et al. (20071 is in agreement 
with the connectivities inferred from our analysis of TRACE and Ha data, 
although the reconnection process may have proceeded in the reverse sense 
compared to that proposed by these authors. 



As stated in Section 3.2.2 it is well-documented that the 1600 A peaks, 
associated with the impulsive phase and B2, constitute signatures of the 
interaction of flare-accelerated particles with the chromospheric footpoints 
of loops connected to the acceleration/injection region. Our results thus in- 
dicate that, during the impulsive phase, particles are preferentially injected 
into LI, while during B2, they are injected into both LI and L2. 



- The time profile of 88345 closely mimics the H4 and H5 light curves. This 
strongly suggests that ^6*345 is optically-thick thermal bremsstrahlung from 
the footpoints of LI. Assuming that an equal amount of 345 GHz emission 
arises from two sources of diameter d (expressed in arcsec), the projections 
of which correspond to H4 and H5, we find that 55*345 leads to a temperature 
increase of up to 8 x 10^ x {20" /d)"^ K. For d = 10" to 20", i.e., comparable 
to H4 and H5 (see Figure [6]), this yields maximum temperature increases 
of 3.2 X 10^ and 8 x 10^ K, respectively. These temperatures are consistent 
with a chromospheric origin for 55*345 and with its similarity to the Ha time 
profiles, the latter being the signature, when seen in emission as during our 
event, of ^lO"' K plasma. 



SOLA: trottet_27Dct2003_FINAL.tex; 28 September 2011; 0:12; p. 22 



Origin of flare time-extended submillimeter emission 



The latter point above and the discussion in Section 3.2.1 suggest that, during 
the gradual phase of the flare, the 230 GHz emission arises from a 1-3 MK source 
of 40" diameter centered near the top of LI, while the 345 GHz emission arises 
from both the same coronal source and from two chromospheric sources of ~20" 
diameter that project onto r4 (H4) and r5 (H5). The convolution of this configu- 
ration of sources with the KOSMA beams (HPBWs of 117 " ±7" at 230 GHz and 
88" ± 5" at 345 GHz, cf. Liithi, Magun, and Miller, |2004[ ) yields source sizes of 
124" at 230 GHz and 107" x 96" at 345 GHz. If we then assume that the emission 
comes from unresolved single sources, after deconvolution we find source sizes 
of 40" at 230 GHz, as expected, and 37" x 61" at 345 GHz, substantially larger 
than the coronal 230 GHz source. Such an increase of source size with frequency 
was previously reported by Liithi, Magun, and Miller (20041, who, during the 
late decay of the gradual phase of the 2001 April 12 flare, measured sizes of 
42" ± 20" and 70" ± 6" at 230 and 345 GHz, respectively. A similar explanation 
as the one proposed here may qualitatively apply to the 2001 April 12 event. 
Indeed, Figure 4 of Liithi, Magun, and Miller ( |2004| ) indicates that the 345 GHz- 
emitting source encompasses both of the two chromospheric footpoints seen in 
hard X-rays, while the 230 GHz source encompasses only the northern one. 



4. Conclusions 

The 2003 October 27 M6.7 flare at ~12:30 UT is one of the few events for 
which a time-extended phase (tens of minutes) following a short (few minutes) 
impulsive phase has been well observed at mm-submm wavelengths. The radio 
data from a few tens of kHz up to 405 GHz shows that the event was confined in 
the chromosphere and low corona, confirming earlier conclusions from the study 
of the active region magnetic field topology. Consistent with the confined nature 
of this event, the 1-345 GHz emission detected during the impulsive phase is 
produced by gyrosynchrotron radiation from high-energy electrons in a dense 
medium. 

The combined analysis of radio total-fiux observations. X-ray imaging and 
spectral data, along with UV/EUV and Ha spatially-resolved observations, has 
allowed us to determine the nature and origin of the radio emission during the 
gradual phase of the flare. The main findings can be summarized as follows: 

- Except for the small impulsive burst B2 where gyrosynchrotron emission 
contributes up to 35 GHz, the radio emission is entirely produced by ther- 
mal bremsstrahlung. 

- Below 230 GHz, the radio emission, which is optically-thin above ~8 GHz, 
is entirely produced in the corona by hot and cool materials at ^--^7-16 MK 
and at ~l-3 MK, respectively. The combination of X-ray data and EUV 
images at 195 and 284 A indicates that the ~l-3 MK plasma fills two large 
coronal loop structures, LI and L2, while the hot X-ray source lies primarily 
within LI, with the centroid of the ~16 MK plasma near the top of the 
brightest portion of the LI loop system. 
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— At 345 GHz, in addition to the coronal emission, there is an optically- 
thick component arising from the lower atmosphere. The similarity between 
the time-evolution of this submm excess and that of the Ha footpoints of 
LI suggests that the excess 345 GHz flux density excess is emitted by 
chromospheric material at a few 10^ K located in the LI footpoints. 



Importantly, these results show that the chromospheric response to the energy 
release during the gradual phase produces negligible emission over a substan- 
tially-extended frequency range, from ^10 to ^230 GHz for the studied flare; 
similar behavior has been observed in other flares, as well. The thermal radio 



emission expected from semi-empirical flare models (cf. Machado et al, 1980 



Mauas, Machado, and Avrett, 1990) is inconsistent with the non-detection of 



chromospheric radio signatures over such a wide frequency range. Indeed, since 
the density-temperature structure considered in these models varies quasi-mono- 
tonically with height, one would expect the radio flux density to increase with 
frequency in the spectral range where the corona is transparent. Thus, observa- 
tions of the radio continuum at cm-submm wavelengths and in the far infrared 
domain can serve as powerful diagnostic tools to investigate both the structure 
and the dynamics of the low solar atmosphere during flares. 
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